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ABSTRACT 



We consider the dissipation by Fermi acceleration of magnetosonic turbulence 
in the Reynolds Layer of the interstellar medium. The scale in the cascade at 
which electron acceleration via stochastic Fermi acceleration (STFA) becomes 
comparable to further cascade of the turbulence defines the inner scale. For 
any magnetic turbulent spectra equal to or shallower than Goldreich-Sridhar this 
turns out to be > 10 12 cm, which is much larger than the shortest length scales 
observed in radio scintillation measurements. While STFA for such spectra then 
contradict models of scintillation which appeal directly to an extended, contin- 
uous tu rbulent cascade, such a se para tion of scales is con s istent with the recent 
work of iBoldyrev fc Gwinnl (120051 ) and lBoldyrev fc Konigll (120051 ) suggesting that 
interstellar scintillation may result from the passage of radio waves through the 
galactic distribution of thin ionized boundary surfaces of HII regions, rather 
than density variations from cascading turbulence. The presence of STFA dissi- 
pation also provides a mechanism for the no n-ionizing heat source observed in the 
Reynolds Layer of the interstellar medium ([Reynolds et al.lll999l ). STFA accom- 
modates the proper heating power, and the input energy is rapidly thermalized 
within the low density Reynolds layer plasma. 



1. Introduction 



Radio scintillation has long been associated with interstellar turbulence (e.g. Minter 
and Spangler 1997)). A major requirement of turbulence-based scintillation models is that 
the inner scale of the cascade is comp arable to the smallest sca l e of the scintillation. An 
altern ative model has been proposed (IBoldyrev fc Gwinnl 120031 . 120051 : IBoldyrev fc Konigl 
20051 ) in which the scintillation instead results from index of refraction variations at the 
photoionized surfaces of a non-Gaussian distribution of either warm ISM regions or HII 
clouds. We find that Stochastic Fermi Acceleration (STFA) is a sufficiently efficient process 
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in the Reynolds Layer that it imposes a cutoff at a scale too large to be consistent with radio 
scintillation, supporting the non-Gaussian cloud model. 



Observations of line ratios by lReynolds et al.l (119991 ) using the Wisconsin H-alpha Map- 
per (WHAM) indicate that the Reynolds layer of the Milky Way interstellar medium (ISM) 
cannot be heated solely by photoionization. Variation with galactic latitude, \z\, of [SII\/Hoc 
and [Nil]/ Ha both are found to be consistent with increases in electron temperature, T e , 
but not in the ionization fractions of eit her species. Additionally, studies of ext ernal galaxies 
NGC891, NGC 4631, and NGC 3079 faevnolds et.alihoool : lotte et.alJ lioOll ) show similar 
evidence of supplemental, non-photoionizing heat sources in the Reynolds layer, the region 
of low density, strongly ionized gas located ^ lkpc away from the galactic midplane. These 
studies supplement the [SI I)/ Ha and [Nil]/ Ha ratios with analysis of [01 1 1] /Hp and 
[OII]/H a , which also are consistent with non-ionizing heating. The heating rate in the 
Milky Way is consistent with being either proportional to electron density, n e , or density 
independent. The heat source is empirically found to have a power input of e = Gin e , where 



G\ ^ 10 erg s , or e = G2, where G2 ^ 10 erg s cm . This allows the supple- 
mental heat source to dominate at high \z\, where n e is low and photoionization heating 
to dominate at low \ z\ w here n e is high. A number of potential mechanisms are presented 



by iReynolds et al.l (119991 ) as supplemental heat sources, including photoelectric grain heat- 
ing, Coulomb collisions with cosmic rays, magnetic field reconnection, and dissipation of 
supperbubble driven magnetic turbulence. 

We consider the role of STFA in the dissipation of interstellar turbulence, looking pri- 
marily for the inner scale of the cascade. ISM turbulence is driven principally by superbub- 
bles, the large blast shells carved out when the members of an OB association reach the 
supernova stage in a short time period. STFA occurs when electrons traveling in along a 
magnetic field line encounter moving compressions of the local magnetic field and are re- 
flected. While STFA drains turbulence primarily into ions for thermal pressure dominated 
plasmas, electrons ar e the primary ener g y recipient when the io n speed is sub-Alfvenic (i.e. 
LaRosa. et al.l (119961 ); iBlackmanl (119991 ); ISelkowitz fc Blackmanl ( 120041 ) ). The latter case is 
relevant for the pre sent study because the Re ynolds layer seems to be a weakly magnetically 
dominated plasma (IMinter fc Spanglerlll996l ). 



Longairlll994l ; ISelkowitz fc Blackmanl 120041 ) 



Fermi 


1949; 


Achterbere 


1981 



2E f 2 

— {v A ±v e v A , 



where E is the initial energy of the reflected particle, c is the speed of light, va is the Alfven 
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speed, and v e is the electron speed. Typical Reynolds layer temperatures are 0.6 < T 4 < 1.2 
where T 4 = T e /10 4 K. This corresponds to an energy in the range 0.5 < E < l.leV. Thus 
5E is significantly below ionization energies, which are tens of eV, and even if the reflection 
rate is high, STFA of protons is a non-ionizing process in the turbulent Reynolds layer of 
the ISM. 

In section 2 we show that STFA provides the correct power and that the energy input 
is quickly thermalized, consistent with the conditions imposed on the heating source. We 
then determine the truncation scale for the turbulent cascade. In section 3, we discuss the 
implications of the inferred STFA dissipation scale for models of interstellar scintillation, 
finding that STFA truncates the cascade at a long length scale, inconsistent with turbulence 
models of radio scintillation. We conclude in section 4. 



2. Stochastic acceleration in the turbulent Reynolds layer 

The power available for electron heating in the MHD turbulent cascade, e-r, can be 
estimated as 



e T 



P 



(2) 



where n e is the electron number density, v a is the loc al Alfven speed, and L is th e outer 
scale of the interstellar turbulenc e inertial range. Fr o m IMinter and Spangler I (119971 ): va = 



2.3 x 10 6 cm s 1 and L = 10 19 cm. IMinter & Spanglerl (119961 ) infer thi s value for L by analysi s 



of emission measure and rotation measure structure functions. From iReynolds et al.l (119991 ): 



n e = 0.28cm 3 at galactic la titude \z 
the hea ting rate called for in 



power (IMac Low, fc Klessen 



26 erg s 1 cm 3 , consistent with 



Revno 



12004: 



lkp c € T = 5 x 10 
ds et al. Jl999h and calcula tions of superbubble injected 



Elmegreen fc Klessen! 12004 ) . 



Selkowitz fc Blackmanl (120041 ) examined STFA in the va 3> c s regime with applications 



to solar flares, and concluded that the post-acceleration spectrum is dominated by the rate of 
escape from the acceleration region. The post-escape spectrum is important in flares because 
the observed x-ray emission is generated via Bremsstrahlung as electrons escape the accel- 
eration region and encounter the dense chromospheric plasma. In the turbulent Reynolds 
layer, the acceleration region is effectively infinite in extent. Instead of observing emission by 
escaped electrons, the observed emission comes from electrons which remain in the Reynolds 
layer. In this case we are concerned with the spectrum of ele ctrons still confined within the 
acceleration region. Regardless, the acceleration rate found in lSelkowitz fc Blackmanl (120041 ) 
remains valid 
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where 

(8E) = Am e v% (4) 
is the average energy per reflection and 



is the rate of reflections, such that 



2L\ L 



(5) 



~dt). 
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and A s is the turbulent inner scale (associated with the parallel component of the magnetic 
field gradient for anisotropic turbulence), Vd is the effective relative speed of the electrons 
and compressions within the plasma, m e is the electron mass, and a is the spectral index of 
the turbulent cascade 

A _ (5B\ a 
L ~ \Bj ■ 

It is assumed throughout our analysis that at length scales comparable to L the strength 
of the magnetic fluctuations, 5B is comparable to that of the mean magnetic field, B. The 
magnetic spectrum of turbulence in the ISM is difficult to measure. Between 10 19 pc and 
0.03pc the spectrum appears to be clos e to K olmogorov (a = 3) (and flatter on larger scales) 
( iHan et al.l (12004] ); iMinter fc Spanglerl (119961 )). There is no direct measure of the magnetic 
spectrum on scales below 0.03pc. 

For the range of magnetic spectra 2 < a < oo the result from our calculations to 
follow, that the cascade truncation scale well exceeds the smallest scintillation scale, will 
not change. Howev er, we first consider a Goldre ich-Sridha r (hereafter GS) turbule nt spec- 
trum with a = 2 (IGoldreich and Sridharl 119971 ) (see also iMatthaeus et al.l (Il998l )). Like 
Matthaeus et al. ( 1998b . GS turbulence is fully anisotropic, and is based on an incompress- 
ible cascade. It involves a more rapid cascade in the direction perpendicular to the local mean 
field than parallel to it. Turbulence becomes less and less compressible on smaller scales, and 
since the Reynolds layer seems to be modestly magnetically dominated (IMinter fc Spangler 
(11996|)), GS is plausible on small enough scales: The GS spectrum more closely arises in com- 
pressible and incomp r essibl e simulations when an initially relative strong field is imposed (e.g. 
Maron fc Goldreich I (hoOlhlCho et al.l J2002ahrVestuto et al.l J2003h Irlaugen fc Brandenburg 
(l2004bl )). i.e. when the initial ratio of thermal to magnetic pressure, (3 < 1 and the velocity 
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fluctuations are of order the initial field. Other values for a should not be ruled out how- 
ever, because the magnetic turbulent spectrum for galactic ISM conditions is not a solved 
theoretical problem. In particular, when an initially weak field is imposed, driven incom- 
pressible simulations show that the magnetic spectral index may be flatter, and flatter than 
the ve locity spectra. Altho ugh it may approach 5/3 for magnetic Prantdl number of order 
unity (IHaugen et al.l (120031)) there may be a trend toward further flatt e ning at larger mag 



netic P randtl number s (ISchekochihin et al.l (120021 ) ; iMaron et al.l (120041 ) ; ISchekochihin et al. 



(120041 ); IHaugen et al.l (12004a ) ). The magnetic spectrum is steeper in the pres ence of kinetic 



helicit y (IMaron fc Blackmanl (120021 ) ) . The compressible driven simulations of IVestuto et al. 



(120031 ) show that the stronger the initial field, the closer the magnetic and velocity spectra 
match, and the flatter the magnetic spectr a the weaker the init ial field. The role of the 
magnetic Prandtl number is hard to assess in lVestuto et al.l (120031 ) since no explicit viscosity 
or resistivity is used. 



Proceeding with a = 2, we note that for this case, ISelkowitz fc Blackmanl (120041 ) find 
that the parallel cascade law is relevant to STFA, and A in eq. (J7j) corresponds to the parallel 
cascade law. When electrons are able to freely stream from one STFA scattering site to the 
next without deflection by pitch angle scattering (the free streaming limit) v a is equal to the 
electron speed v e and, 



\dt) f 



L 



-m P v A v P 



L 



-1/2 



In order to determine the truncation scale of the cascade, we impose the balance condition 
that the turbulent power, er is equal to the STFA acceleration rate. Setting the STFA rate 
equal to the total power selects a single value of A/, which is an effective upper limit on the 
truncation scale. We use the subscript / to denote the free streaming limit. 

In the limit of very strong pitch angle scattering, where the scattering length scale 
\ p <C A/, A p can be considered the electron mean free path. Electrons are effectively trapped 
in regions smaller than A/, executing a random walk with small drift speed. They encounter 
magnetic compressions only as the compressions stream past, significantly reducing the ac- 
celeration rate below that of Eq [HJ The power balance condition between STFA and the 
cascade power is not met, and the cascade will proceed to scales shorter than A/. The cas- 
cade cannot continue farther than a length scale comparable to X p . At this point, the system 
passes back into a free-streaming regime, with a substantially shorter A than is required to 
meet the power balance condition. The turbulence then drains rapidly. 

From the above two paragraphs, we conclude that the cascade would be truncated at a 
scale no smaller than the minimum of A p and A/. To assess the minimum scale of dissipation 



- 6- 



for the Reynolds layer, we ass ume that the pi tch angle scattering is dominated by electron- 
electron Coulomb scattering. ISpitzer I (119671 ) finds the Coulomb self-collision time, t c , to 
be 



n e in A 

where T e is the electron temperature, and In A ^ 25 is determined by the effective long range 
cutoff of the Coulomb force in a plasma. The electron speed is given by v e = ^/2kT/m e , 
and thus 



\<n Vptr 



2&T 0.266 T f 



3/2 



n„ 



In A 



2 x 10 id T, 



13^2 
4 i 



(10) 



where T 4 = T e /10 4 .For the range of observed temperatures in iReynolds et al.l (119991 ). 0.6 < 
T4 < 1.2, the electron mean free path falls in the range 8 x 10 12 cm A p < 3 x 10 13 cm. Taking 
the lower limit of this range fixes the lower bound of the cascade truncation scale. 

In the free streaming limit we can find Xf by setting = ;/ 



dE \ 
e ^ dt >fi 



A / =f^Vf^VL = 6xlO"ir l . 
\m p J \v A J 

For the observed temperature range of 0.6 < T4 < 1.2, 10 15 cm> A/ > 5 x 10 14 cm 
an upper bound on the cascade truncation scale. 

It should be noted that the free streaming approximation is invalid in the Reynolds 
layer as A//A p = 100 at T 4 = 6. There are, on average 100 pitch angle scattering events 
per electron per encounter with a compression. Given the low fraction of encounters which 
result in a reflection, F, electrons cannot stream freely from one reflection site to the next. 
The quantity F is determined by the pitch angle condition for reflection 



2n SB 

cos 9 min < — , (12) 

£5 

and the turbulent cascade law (e.g. eq. ((Tj) for GS turbulence). All electrons with 8 > 9 min 
reflect. Assuming pitch angle isotropy, we estimate 



F = i- f 2vrsin(^=i(A / /L) 1/4 , (13) 

where we have integrated over all angles greater than Q m in to find the fraction of phase space 
which satisfies the pitch angle condition for reflection, and we have used (J7|) for the last 
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equality. Electrons which have too large a component of their momentum in the direction 
parallel to the field are not stopped by the compressions, and pass through them. The 
acceleration rate must be retarded significantly. We have bounded the turbulent dissipation 
scale 8 x 10 12 cm< A s < 10 15 cm for GS turbulence for a Goldreich-Sridhar, a = 2, cascade. 

In order to produce a dissipation scale consistent with the smallest scales observed 
in scintillation measurements, the cascade must truncate at the inn er scale predi c ted b y 
scintillation models. Stated value s of this scale vary: 3.5 x 10 6 cm in iMoran. et.al.l (119901 ). 
3 x 10 7 cm in iMolnar. et.al.1 (119951 ). and 3 x 10 10 cm in [Rickety (119901 ). Rewriting eq [TT1 for 
an arbitrary value of a, we have 



A, 



2m f 



l-l/a 



VA 



l-l/a 



(14) 



In order to produce a value A/ = 10 10 cm, the cascade must have a very steep spectrum, with 
a < 1.3. The other predictions require even steeper spectra. Such spectra are inconsistent 
with simulations of interstellar turbulence, as well as the inferred spectral indice s of scintil- 
lation models, which place the lo wer limit near a Kolmogorov spectrum, a = 3 (Han et al. 
20041 ; iMinter and Spanrier Ill997l ). 



Assuming the Reynolds layer is weakly magnetically dominated (IMinter fc Spangler 



19961 ) implies t hat we have studied di s sipati on of fast magnetosonic mode turbulence. On 
the other hand, iLithwick fc Goldreichl (120011 ) argue that for thermally dominated plasmas, 
(3 > 1, the turbulent density fluctuations are dominated by the slow and entropy mode, and 
that the fast mode is decoupled from the slow, entropy, and Alfven modes. Of the slow and 
entropy modes, only the slow mode is important for STFA. While the fast mode could be 
dominant for STFA in the Reynolds layer, the acceleration by slow modes does achieve the 
same qualitative result, albeit at a shorter dissipation scale. In the Reynolds layer, where 
(3 ^ 0.1, the slow mode velocity is given roughly by the sound speed c s ^ 0.2^. If both 
modes are present with equal energy density, the fast mode is more efficient for STFA than 
the slow mode. To apply eq. to STFA of slow modes in a plasma with (3 < 1 requires 
the replacement of va with c s , such that 



dE\ 2 2 /A s y (1 - 1/a) 

—m e c s v d [ — . (15) 



dt J s L \L 

For free streaming electrons, there is no change in Vd = v e and the dissipation scale is only 
decreased by a factor of (c s /va) 2 = 0.04. For trapped electrons, the mean free path continues 
to place a lower bound on the dissipation. Despite the moderate decrease in X s in a slow 
mode dominated cascade, it remains well above the scintillation length scale. 



- 8- 



In addition to providing the correct power input and a physically plausible dissipative 
scale, STFA must sup ply energy in a manne r consi stent with heating, as opposed to non- 



thermal acceleration. ISelkowitz &: Blackmanl (12004 ) demonstrated that, in the limit of no 



self-interaction which is appropriate for solar flare applications above a few keV, the electron 
energy spectrum within the acceleration region is shaped by two competing effects: a bulk 
shifting to higher energy, and a diffusive spreading of the spectrum. The resulting spectrum 
can be quasi-thermal, even when electrons are not truly sharing energy. However, in the case 
of ISM heating, where Coulomb self-scattering of electrons is presumed to be the dominant 
source of pitch angle scattering, STFA is accompanied by rapid thermalization of the electron 
population. Consider the acceleration time, t S tfa 

Tstfa = yjy- = 3 x 10 12 T 4 s, (16) 

\ dt ) s 



and the thermalization time, r eq , ( jSpitzer Ill967l ) 



= 6 x 10 5 T 4 3/2 s, (17) 

where we have taken the thermalization time to be the two species equilibration time, with 
both species being electrons at a single temperature, T 4 . When r eq < tstfa, as is the case 
for the Reynolds layer, then the energy input via STFA is shared among electrons rapidly; 
the energy spectrum is thermal. The electron-proton thermalization time is greater than r eq 
by a factor of the mass proton to electron ratio m p /m e = 1836, which is nevertheless shorter 
than the acceleration time. The protons (and heavy ions) remain in equilibrium with the 
electrons as well. 



Implications of the turbulent cutoff scale for models of interstellar 

scintillation 



The dom inant power source fo r inter s tellar turbulence is lik e ly sup ernovae and super- 
bubble shells (jElmegreen &: Klessenl ( 120041 ) ; lMac Low, fc Klessen I ( 120041 ) ). Superbubbles are 
powerful and frequent enough to periodically pass through, and thus reseed turbulence in, the 
entire galactic disk. If the efficiency for converting a single supernova's mechanical luminos- 
ity of 10 51 erg to turbulent energy is t]sn = 0.1, the rate of supernova events in the Milky Way 
is Rsn = 50yr _1 , then the total turbulent driving from supernovae and superbubbles is esti- 
mate d to be 3 x 10~ 26 erg s -1 (reviews by (IMac Low, fc Klessen 1120041 ; lElmegreen fc Klessen 
20041 ) and references within). This is consistent with t he estimated damping rate from eq 
(1) as well as the required supplemental heating rate ((Reynolds et al. 1999 ). 
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Although we have considered a GS cascade law because it is inherently anisotropic, and 
incorporates the magnetic field into the formalism directly, we emphasize again that the GS 
cascade, with a = 2 magnetic spectral i ndex li es at one end of a r ange of possible cascade 



power laws (e.g. iGoldreich and Sridharl (119971 ); I Vestuto et al.l ( 120031 )). As discussed, varying 



a over a wide range does not substantially alter the results of our analysis. 

Likewise, it is n ot certain which magn e toson ic wave mode, fast or slow, is dominant in 
the Reynolds layer. iLithwick fc Goldreichl (120011 ) argues for slow modes to dominate in a 
high (3 plasma. Fast mode turbulence may play a stronger role in (5 ^ 0.1 diffuse Reynolds 
layer plasma. Because the fast mode is a more efficient STFA accelerator, due primarily to 
its higher phase speed, we emphasized STFA on fast modes, but it was shown in section 
2 that the major qualitative result of our study holds for either wave mode: STFA is a 
sufficient source of heating, and the dissipation scale of MHD turbulence is longer than the 
short scales observed in radio scintillation measurements. 

The correspon dence between t h e turb ulent driving power and the inferred supplemen- 
tal heating rate of iReynolds et al.l (119991 ) is suggestive of a connection. This connection 
would be more strongly verifiable if the turbulent dissipation scale could be linked t o di- 
rect observations of the Rey nolds layer. This has not yet been done. ISpangler I (|199ll ) and 
Minter and Spangler I (119971 ) considered that the radio scintillation observations (e.g. Arm- 
strong, Rickett, Spangler 1995) may be consistent with scattering by localized turbulent 
density fluctuations across the inertial range of Kolmogorov turbulence, but they further 
suggest that the "fluctiferous" medium is either HII region envelopes, or the higher density 
portions of the warm ionized medium (WIM). Although possible turbulent dissipation scales 
(like the Larmor radius) of the scintillating region correspond well with the smallest scales 
of radio scintillation, no scintil lation measurements are ana lyzed in lines of sight through 
the low density Reynolds layer. iMinter and Spangler I (119971 ) instead assume that either the 
Larmor radius or the ion inertial length is the key dissipative scale in the Reynolds layer, 
just as in denser regions of the ISM, and use this assumption to calculate dissipation and 
heating rates. 

We find that the STFA turbulent dissipation scale is instead likely far larger (A s > 
10 13 cm) in the Reynolds layer, for GS turbulence. The STFA dissipation scale does drop 
to the observed scintillation scale in the case of a very steep turbulent spectrum at small 
sc ales, with index a < 1 . 3. In terpretations of the scintillation as turbulence driven, such 



as 



Minter and Spangler I (119971 ). indicate a spectrum no steeper than Kolmogorov, 



= 3. 

Although STFA is seemingly at odds with radio scintillation measurements, it may not be if 
the source of the scintillation is not a continuous turbulent cascade as is commonly assumed. 
In particular, recent studies have found that interstellar radio scintillation is in fact consistent 
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with a non-Gaussia n distr ibution of stellar ionization shells iBoldyrev fc Gwinnl (120031 . 120051 ); 
Boldyrev &: Konigll (120051 ). Local index of refraction variations at the thin photoionized 
surfaces of molecular clouds are distributed randomly throughout the interstellar medium. 
Radio waves are deflected through small angles at these fronts, and adjacent rays can be 
deflected through differing paths. A detailed model employing Levy statistics provides an 
excellent fit to the temporal spreading of pulses from dis tant pulsars, which tradit ional 
cascade scintillation mode ls have had difficulty explaining. IBoldyrev fc Gwinnl (120051 ) and 
Boldyrev fc Konigll (120051 ) conclude that the density fluctuations are more likely to result 
from randomly distributed thin ionization shells in the ISM than from global turbulence. 
As a result, the scintillation would no longer provide a constraint on the inner scale of the 
turbulent cascade. Turbulence driven on very large (> 10 pc) scales could damp on scales 
larger than the smallest scintillation scales and dissipation of the turbulence by STFA would 
not conflict with scintillation measurements. 

The low er bound of the dissipa t ion sc ale found above violates an assumption used in the 
approach of iMinter and Spangler I (119971 ). Given that the measured supplemental heating 
rate in the low density ISM can be well explained by turbulent dissipation, one may infer 
that turbulence is driven and dissipated throughout the ISM. The density is higher and 
temperature lower in the lower \z\, and thus X p is shorter, in the portions of the WIM where 
the scintillation may occur. STFA is then even more closely tied to the lower bound for 
A s set by Coulomb scattering at low \z\. If STFA is important, it would then be unlikely 
that turbulence input by superbubbles cascades to sufficiently small scales in any region of 
the interstellar medium to account f or the scint il lation observations. This g ives credence to 
the argument of IBoldyrev fc Gwinnl (120031 . 120051 ); IBoldyrev fc Konigll (120051 ). which does not 
rely on cascading turbulence but concludes that the statistical spatial distributions of the 
ionized boundaries of molecular clouds better explains scintillation than do models based on 
turbulence dominated density fluctuations. 



4. Conclusion 



Dissipation of turbulence almost certainly plays a role in the heating of the interstellar 
medium. The turbulent energy supply from superno vae is sufficient to provide the supple 



(2000 


); 


Otte et.al. 


(2001) 



the Reynolds layer of the Milky Way, truncating the turbulent cascade at length scales no 
shorter than 8 x 10 13 cm for a cascade no steeper than Goldreich-Sridhar. This truncation 
scale is too large to be consistent with turbulence-based models of interstellar scintillation. 



Instead, models like those of iBoldyrev fe Gwinnl (120031 . 120051 ); 
which do not rely on the cascade, are required. 



Boldvrev fe Konigll (120051 ). 
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